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Abstract 

We apply a new bosonization technique to relativistic field theories of 
femiions whose partition function is dominated by bosonic compos- 
ites, and derive the effective action for these bosons. The derivation 
respects all symmetries, including gauge invariance, with the excep- 
tion of Euclidean invariance which must be checked a posteriori. We 
use a lattice regularization which should make applications to gauge 
theories easier. We test the method on a fermion field theory with 
quartic interaction in the limit of large number of flavours Nf, and 
show that it reproduces the exact results in the bosonic sector, namely 
condensation of a composite boson with the right mass which breaks 
the discrete chiral invariance of the model. Moreover we determine 
the structure function of the condensed composite, whose spatial part 
turns out to be identical to that of the Cooper pairs of the BCS model 
of superconductivity. 
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1 Introduction 



There are many theories involving fermions whose low energy excitations 
are dominated by bosonic modes. Historically the first case of great general 
relevance is superconductivity in metals [T], which was the starting point 
of innumerable applications and theoretical developments in several fields. 
In the context of relativistic field theories important examples are vector 
dominance in strong-electromagnetic interactions [2] and dominance of chi- 
ral mesons in QCD. In these cases the dominant composite bosons have 
fermion number zero, but in QCD at high density so called colour supercon- 
ductivity [3j is expected to occur, and in this case the dominant composite 
bosons have fermion number 2 like Cooper pairs in superconducting metals. 
Obviously there is an immense literature about these subjects, but since we 
address the problem of composite boson dominance in its generality without 
actual applications to any of the above examples, we quote only historical 
works. We note however, that in spite of a huge literature about bosoniza- 
tion all these systems are, so far, mostly treated phenomenologically. 

Nonrelativistic methods can hardly be extended to relativistic field the- 
ories with few exceptions, one being ref. 0] which has some features in 
common with the relativistic papers quoted below. We therefore only men- 
tion Bogoliubov's work on superconductivity for his historical importance 
as the first organic approach to bosonization |,5j. We do not discuss fermion 
systems in 1+1 dimensions either. For them a wealth of exact results exists 
which are however dimension specific and can be found in textbooks [6]. 

For relativistic theories in higher dimensions the problem of bosoniza- 
tion has essentially two features: how to introduce the fields associated to 
the dominant composite bosons and how to handle the resulting effective ac- 
tion. There are basically two approaches. In the first one higher dimensional 
terms (quartic in the fermion fields) if absent, are introduced by hands in 
the action and the bosonic fields are generated by the Hubbard-Stratonovich 
transformation [7]. Sometimes, when the added terms are irrelevant in the 
renormalization-group terminology, they are introduced to stabilize the eval- 
uation, for instance in numerical simulations, of the fermion determinant [8], 
but, more often, the investigations are analytical studies and are substan- 
tially restricted to a nonrenormalizable framework In one noticeable 
exception [9] the Hubbard-Stratonovich transformation is scale-dependent 
("re-bosonization") leading to a functional renormalization-group analysis. 
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In the second approach the boson fields are Kalb-Ramond fields introduced 
by a change of fermion variables in the action [TOj , a technique also adopted 
in the quoted work [4j about non relativistic many-body systems. No general 
procedure however has been developed to handle the fermionic determinant 
appearing in the effective action, in particular for gauge theories. Moreover 
the above approaches are essentially restricted to composite bosons with 
fermion number zero. 

Recently a method of bosonization was developed in the framework of 
many-body theories by which we can treat both charged and neutral com- 
posites [11]. The starting point in this method is the partition function in 
operator form, namely the trace of the transfer matrix in the Fock space 
of the fermions. The physical assumption of boson dominance is then im- 
plemented by restricting the trace to fermion composites. This requires an 
approximation of a projection operator in the subspace of the composites, 
the approximation being the better, the higher the number of fermion states 
in the composites. The approximate projection operator is constructed in 
terms of coherent states of composites, and evaluation of the trace, which is 
done exactly, generates a bosonic action in terms of the holomorphic vari- 
ables appearing in the coherent states. 

Our approach shares two features with a variational method: The reduc- 
tion of the starting full space, here the fermionic Fock space, to a subspace, 
here that of the composites, and a variational procedure to determine the 
structure functions. The utility of variational methods and bosonization 
has been widely appreciated in the theory of many-body systems. But their 
potentiality has also been considered in the framework of relativistic field 
theories, in particular gauge theories, for for example by R. Feynman [12j 
who, however, was sceptical about its practical applicability, and recently in 
connection with QCD at high baryon density [13j . 

The particular bosonization method we investigate has been completely 
developed for nonrelativistic many-body systems and checked on the BCS 
model of superconductivity and the pairing model of finite systems like 
atomic nuclei and ultra-small superconducting grains [11]. The properties of 
these systems are exactly reproduced. But in the nonrelativistic derivation 
of the effective action, higher powers of the temporal spacing, which would 
not contribute in the continuum limit, were neglected. This obviously can- 
not be done in relativistic theories because of ultraviolet divergences and in 
the present paper we perform a different derivation. 

We adopt a lattice regularization for two reasons. The first one is that it 
allows an unambiguous definition of composites. The second one is related 
to gauge theories. We expect that in these theories the action of effective 
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bosons will involve vacuum expectation values of invariant functions of gauge 
fields which cannot be evaluated in the present framework. A lattice for- 
mulation should allow us to extract such expectation values from numerical 
simulations. 

There is a price to be paid for such advantages related to the well known 
difficulty with chiral invariance on a lattice, which can only in part be over- 
come by using Kogut-Susskind fermions. However our method can, at least 
in principle, be used with any regularization for which a transfer matrix has 
been explicitly constructed. 

A major difference with respect to nonrelativistic theories is related to 
Euclidean invariance. The formalism of the transfer matrix does not treat 
time and space in a symmetric way, and therefore Euclidean invariance of 
the bosonic action must be checked a posteriori. All other symmetries are 
instead respected in our procedure. 

We tested the validity of our method on a model with 4-fermion interac- 
tion [T3] in 3+1 dimensions: we exactly reproduce all the known results in 
the boson sector, namely condensation of a composite boson which breaks 
the discrete chiral invariance of the model and its mass. In addition we de- 
termine its structure function, whose spatial part is identical with that of the 
Cooper pairs of the BCS model of superconductivity. Our approach in this 
case has some points in common with a variational calculation performed 
for the model in 1+1 dimensions [15j . 

Our solution of the model is more complicated than the standard one, 
but it is not done looking for a greater simplicity, but as a check with special 
attention to Euclidean invariance. In any case it gives explicitly the structure 
function of the composite, which to our knowledge was not known. 

In this work we confine ourselves to fermion systems at zero fermion 
density. An extension to non-vanishing fermion density will be presented 
elsewhere [16] . Its application to the 4-fermion interaction model reproduces 
the properties of the fermion sector. 

The paper is organized in the following way. In Section II we report, for 
the convenience of the reader, the general formalism. In Section III we derive 
an efi^ective action for an arbitrary relativistic theory, and discuss the saddle 
point approximation. In Section IV we present an alternative derivation of 
the effective action with the relative saddle point approximation which turns 
out to be identical to the previous one. In Section V we apply our method 
to the model with 4-fermion interaction and in Section VI we summarize 
our results. Most technical details are relegated in Appendices. 
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2 General formalism 



For reader's convenience we report the general formalism of bosonization 
based on coherent states of fermion composites. By using discrete time, the 
operator form of the partition function is 

ZF = Tv^ll%,t+i. (1) 
t 

The trace is over the Fock space of the fermions and ^,t+i is the Euclidean 
transfer matrix which maps the Fock space at time t into that at time t + 1. 
t runs in a range which depends on the temperature. 

Under the assumption of boson dominance we can restrict the trace to 
fermion bosonic composites. The restricted partition function can be written 

Zc = Tr^l{rTt,t+i, (2) 

t 

where P is a projection operator in the subspacc of the composites. Because 
the formalism treats asymmetrically time and space, it is convenient to use 
the following notation: we shall use boldface letters, as x, to denote spatial 

coordinates, and italic letters to denote space-time coordinates: x = (t, x). 

To construct the projector V we first introduce the composite creation 
operators 

K,K = <^HkV^ = E ^MK)ijv} , (3) 
ij 

where x represents the spatial coordinate of the composite and K its quan- 
tum numbers, among which can be radial excitations, orbital angular mo- 
mentum, spin, flavour, etc. The are the composite structure func- 
tions {wave functions) and uj and vj are, respectively, creation operators of 
fermions and antifermions in state i, obeying canonical anti-commutation 
relations, 

{ul,Uj} = = Sij , {ui, Uj} = {vi, Vj} = {ui, Vj} = {nj, %} = . (4) 

Since the fermion creation operators are nilpotent, the composite cre- 
ation operators can be classified according to their index of nilpotency, 
which is the highest integer exponent such that 

{^^Y^^O. (5) 
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The composite structure functions ^:x_k o-^^ to be determined variationally 
in order to maximize the saturation of the partition function Zq ■ 
It is useful to introduce the operator doublet 



and the orthogonal projectors 





,(+) 



u 



in such a way that 



(6) 

(7) 
(8) 

(9) 



By analogy with canonical bosonic systems, we build coherent states of 
composites 



10 = exp 



|o>, 



(10) 



where the ^xii''s are holomorphic variables and |0) is the fermion vacuum. 
We call these states coherent because they have the form of coherent states 
of elementary bosons, sharing with these states the property of a fixed phase 
relation among the components with different number of bosons. But the 
basic property of coherent states cannot be fulfilled. Indeed 



(11) 



because composite creation-destruction operators don't satisfy canonical 
commutation rules. Again in analogy with canonical bosonic systems, we 
define the operator 



V 



where 



d^dC 

2TTi 



d^dC 
2m 



n 

■>s.,K 



1 



10 (ei 



27rz 



(12) 



(13) 



which is neither a projector nor the identity in the subspace of the compos- 
ites, but, as shown in Appendix [Bl it approximately becomes a projector 
onto the space of composites when the latter ones have a large index of 
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nilpotency. Indeed since the composites do not obey canonical commuta- 
tion relations, their properties can be very different from those of canonical 
bosonic coherent states. However, if the index of nilpotency of the com- 
posites is large enough, the composite system resembles a canonical bosonic 
system, and all the properties of canonical boson coherent states will ap- 
proximately hold for the composite coherent states. 
The scalar product of coherent states is 

(e|0=det+(J + i?i?'t) (14) 

where 

B^ = ^.<^^ = Y,^^K<K (15) 

and for any matrix A we define 

det±A := det(P(j^^A) . (16) 

The projection operator Pq^^ appears when matrices act on half the fermion 
field. Notice that the entries of all the matrices do not include time. X is 
the identity in the space of these matrices. By a little abuse of notation we 
will write "1" instead of I when there should be no ambiguity in the inter- 
pretation. Similarly we will replace by "1" the identity in various subspaces 
like color, taste and so on. 



3 First form of the effective boson action 

According to |18j , the fermion part of the transfer matrix, which is a function 
of all elementary bosonic (also gauge) fields, that we represent by a, is 
factorized as 

Tt^t+i = T [at, at+i] = ft [at] f [at+i] = f} f,^, . (17) 

We use the subscript t to denote the dependence of any matrix on the 
particular configuration of bosonic fields at, and 

f = exp[-u^M u - v^'M'^v] exp[vN u] (18) 

aside from a possible extra factor which is a function only of the bosonic 
fields and can be therefore included in the bosonic contribution to the parti- 
tion function. The form of the matrices M and N (the superscript T means 
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transposed) depends on the nature of the interactions and the regulariza- 
tion adopted for the fermions. What follows does not depend on their ex- 
plicit expressions which are reported in the Appendix [C] for Kogut-Susskind 
fermions in the flavour basis (we do not know any suitable expression in the 
spin-diagonal basis). 

At finite time-strip of length Lq, with periodic boundary conditions for 
the bosonic fields, the partition function ([T]) is 

Zp = Tr^f^f,fl...fl_,f, (19) 

while its restriction to the composites ([2]) is 

Zc = TT^Vf^f,Vfl---Vfl_^fQ 



Lo-l 

n 

t=0 



2TTi 



m-. 



^ (edT/T;+iiem) (20) 



where we have introduced a copy of the Fock space of the composites at 
each time slice. Explicitly 

16) = exp Ck {t, x) [at}j |0) . (21) 

We remark that the structure functions ^ do not depend explicitly on time, 
but as they are defined in presence of an external bosonic field configuration, 
time will enter as a label of the bosonic fields. 

We introduce also copies of the matrix B, defined in (jlSp at each time 
slice 

^l = (^^^)^ = E^/^(*'^)^xKN• (22) 

Now setting M = 0, a restriction which will be eliminated in the next 
Section, we evaluate the matrix elements 

T;+i|6+i) = (0|e^'^'"e"^^* ^^e^^*+i"e"^^*+i«^ |0) . (23) 
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With the help of the formulae collected in Appendix [A] we find 



da* dad (3* d (3' 
{af3\ap) 



{0\e 



a(3){a(3\e 



vNt+iu bI^-^v'^ 



|0) 



det+(l + i^iViiV^^i) det+(l + BtNl 



X / [da* dadl3* dl3]e 



-a'a-lS'lS-fS* 



Bf, ,a*-aBT- 



1+JV*S 



det+(l + 5;+iA^t^i) det+(l + BtNj 
X det_ ^1 + b}^-^ ^ 
det+ A,t+i 



l + Nt+iBl,l + BtN} \ 



where 



Dt,t+i = (l + BtN}^ (l + iVt+iSj+i) + 5*4+1 
Therefore by substitution in ()20|) . using equation ()14p 



n 



27ri 



exp [-5eff(r,0] 



where 



n 



2m 



5, 



eff 



n 

J^tr 



27ri 



In 1 + St 5/ -lnA,m 



(24) 
(25) 

(26) 

(27) 
(28) 

(29) 



In analogy to p6p we used the definition valid for any matrix A 

tr±A :=tr(P(S^^A). 

and we have replaced (t,x) by x to shorten the notation. 

It is remarkable that the effective action for the composites Seff has been 
evaluated exactly, so that the only approximations in the partition function 
are the physical assumption of boson dominance and the form of the projector 
over the subspace of the composites. 

Remember that the entries of these matrices do not not include time, 
which appears only as a label. For instance the matrix D satisfies the 
equation 



Dt 



H,t+i=Dl^^^,, (30) 
where the operation of Hermitean conjugation does not affect time. 
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3.1 The saddle point equations 



When the index of nilpotency is large the partition function is dominated by 
the minimum of the action. Its variation with respect to the matrix elements 
of B and B^ gives 



B 



t 1 
* 1 + BtBl 



NUl + Nt+iBl^,)+Bl^, 



1 



t+i 



-B. 



l + Bt^,Bl^, 



t+i 



D 



t,t+i 



1 + BtNl Nt+i + Bt 



(31) 
(32) 



As a consequence of (|30j) . these equations are not independent from each 
other, because the second can be obtained from the first one by exchanging 
t with t + 1 after Hermitean conjugation. So it is sufficient to study only 
one of them. 

In the sequel we consider the case in which the effective action has a 
minimum with respect to the bosonic fields, coupled with the fermions, at 
a constant value. Of course, this assumption cannot hold for gauge fields. 
We are then able to determine the solutions which provide the composite 
structure functions at the semiclassical level. 

If the bosonic fields appearing in Nt are constant at the saddle point, 
there are constant nontrivial solutions Bt = B for the second order equation 



Setting 



for ^ / we find 



BN^B - BN^N -N = 0. 



B = NA 



N (a N'^NA - A N^N - 1 1 = . 



(33) 



(34) 



(35) 



If we separate A into its Hermitean and anti-Hermitean parts, we see from 
this equation that A commutes with N^N. The solution 



A, 



where 



2H 



ij2 ^ _jvtiv. 
4 



(36) 
(37) 



shows that A is Hermitean. 
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Remark that 

-1 



(38) 

A+ + A_ = 1 . (39) 
The effective action at the saddle points 

5± = - ^ tr+ In {l + iVA±iVt) (40) 
t 

= - 5^ tr_ In (^A\N^N^ (41) 
t 

\±2 



^ tr _ In y/l + H^ + j (43) 



takes opposite vahics, the minimum being for the upper sign. The ground 
state is a condensate of composites whose structure function in polar repre- 
sentation is 

B = ^NH-^f{H), (44) 

the polar radius being 



f{H) = Vl + H^ + H. (45) 

In such composites the occupancy of high momentum fermion states is larger 
than that of low momentum states. Such a structure of condensed bosons 
is quite different from the structure of (even virtually) bound pairs. But it 
can be set in a more natural form by the unitary transformation 

{v')^ = u, («')^ = V . (46) 

It defines a new vacuum 

u'\0') = v'\0') = (47) 
related to the original vacuum according to 

\0') =l[ulvl\Q) . (48) 

i 

The new vacuum is the trivial solution of the saddle point equations. It is 
the completely filled state and it is then physically equivalent to the original 
vacuum. 

In the next Section we will derive a new form of the effective action using 
the new creation-annihilation operators. 
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4 Second form of the effective boson action 



We perform the evaluation of the effective action with the new operators. We 
first remark that the transformation (j46p interchanges the role of the projec- 
tors Pq^^ ■ Since the eigenstates of these projectors correspond to fermions 
propagating forward, respectively backward in time, this transformation is 
related to time-reversal. Under its action 



(49) 



In the general case this transformation induces the replacements 

N ^ N'' 

M ^ -M (50) 

and a change in the purely bosonic contribution to the action. At M = 
it interchanges in form T with whenever N is Hermitean. Now (|20p 
becomes 



2c 



Tr^ VfQflvf-^---V f^g_i f(J 



Lo-l 

n 

t=0 



2m 



(51) 



where we have used the same notation as before for the coherent states but 
now 

^|0'). (52) 
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To construct the effective action we have to evaluate the matrix elements 



(O'l 



t 



exp [^t+iv'^^u') n 



10') 



(53) 



Now we notice that for an elementary boson the basis of coherent states 

exp(6 6)(6^)^|0), exp(66t)|o), (54) 

are equivalent to each other in the subspace of states with 0,1, ... ,Q bosons. 
By analogy instead of the states exjp{^v'^^u') x ni(''^D^(^i)^|0') '^^ can use 
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the basis exp {v')^) |0'). They are not exactly equivalent because the 

composites do not satisfy canonical commutation relations, but they should 
be equivalent within our approximations, and we will have a crosscheck 
of this assumption. In the following we drop the prime on all creation- 
destruction operators. 

First we evaluate the action of on coherent states 

f t|0 = exp['u^iVf exp[-{i^Mt u - v^^M*!)] exp[u^ B'' v^O) . (55) 

By repeated application of the identity, valid for arbitrary matrices A, B 

exp[ii1'Su]exp[nt yl?)t]|o) = e^Y>[v) Av^]\id) (56) 

we get 



exp < u 



V 



}|0) 



(57) 



Then the matrix elements of the transfer matrix are 

(61 f, = det+ {1 + [Nt + eM-Mt)Bt exp(-Mt)] 

X [nI, + eM-Ml,)Bl, exp(-M/+i)]} (58) 

and the new effective action reads 



cS 



1 + BtB^ 



t 

x(<i + e<^<,e<.)]} 



ln[l+(7V, + e^*i?ie^') 



in which we have changed the sign of the matrix M according to (jSOp . 
Notice that while the partition function must not change by application of 
the unitary transformation (j46|) . the effective action need not remain the 
same. For our application to the model with 4-fermion interaction we are 
interested in the case of Mt = 0. In the following we shall restrict ourselves 
to this case, then 



S',s = Y.tT+ \ln(l + BtBl)- In D[ 



t,t+i 



with 



D',^,^, = l + {Nt + Bt) (nI^ + bI, 



(59) 
(60) 
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Let us investigate what is the relation between this effective action and the 
one we got in the previous section. Its expression (j28p can be rewritten as 

t 

-InjSt [{b^^ + n}) ((i?l+,)-i+iVi+i) +i] bI^]] 

= ^tr_ {ln{l + B^HBl)-') 
t 

-In [l+(i?-i + Art) ((i?t^,)-i+iV,+i)]} 
which coincides with under the change 



Bt 



Bt 



B}Bt 



(61) 



for every t. Under this change the saddle-point equation for time indepen- 
dent solutions Bt = B = NA becomes 



N^NA^ + N^NA^ - 1 1 = {) . 



(62) 



This is the same as the equation we had in the previous section under the 
change A — > —A^ therefore we have the time independent solutions 



B± 



N 



-NA± = ^ — {Vl + H^±H). 



(63) 



We see that in the saddle point approximation <I> is Hermitian when is 
Hermitian. At the saddle point 



^tr_ In (l + A^N^N 



A, 



tr_ In (Vl + H^ + H 



T2 



takes its minimum for the lower sign, and 

S+ = s'_. 



(64) 
(65) 
(66) 

(67) 



At the minimum the values of the two actions coincide, and therefore the 
corresponding partition functions are equal at leading order. Let us notice 
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that the transformation ([6T]) . because of the identity ([38|) . interchanges 
with —A-. 

The trivial and nontrivial solutions are strikingly similar to the corre- 
sponding solutions of the BCS model of superconductivity. The comparison 
is best done using a polar representation for the structure functions which 
separates the unitary factor which depends on the intrinsic degrees of free- 
dom, from the polar factor which is a function of the spatial coordinates. 
Now the unitary factors are necessarily different, because different are the 
intrinsic degrees of freedom, but the polar factors and therefore the spatial 
wave functions turn out to be identical. Indeed in the case of superconduc- 
tivity the trivial solution is given by = oo for k < kp, = 0, for k^ > kp, 
k being the momentum and kp the Fermi momentum of the electrons while 
the structure function of the condensed Cooper pairs is given by [TT] 



The Pauli matrix cj2 which couples the spins to zero is replaced in a rela- 
tivistic theory by the spin-taste structure matrix A^(2 H)~^ and the electron 
kinetic energy E measured with respect to the chemical potential divided 
by the energy gap is replaced by the relativistic energy H. In the present 
case the chemical potential disappears, because the composites are neutral, 
and to investigate the effect of a chemical potential one should introduce 
states with non vanishing fermion number. What makes the comparison 
with the BCS model somewhat involved is the absence of a Fermi energy in 
relativistic theories with zero chemical potential. Indeed the fully occupied 
state in the old operators of the relativistic theory corresponds to the state 
where all single particle states are occupied below the Fermi surface in the 
BCS model. 

5 Application to a model with 4-fermion interac- 
tion 

We apply our formalism to the field theory with quartic interaction in 3+1 
dimensions regularized on a lattice with Kogut-Susskind fermions. For each 
of the four Kogut-Susskind tastes we take Nj degenerated flavours. Hence, 
the continuum limit will describe a theory with 4Nf flavours. In the flavour 
basis the action reads 



'5 = E ' E '^(^) ["^ ^ ® ^ + Qky V'(y) + 2 ^ E (69) 




(68) 




X 



y 



X 
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where m is the mass parameter, the couphng constant, ip the fermion 
fields and Q the hopping matrix: 

q = J2j^01 [p^-)vW + pWv(-)] . (70) 

The matrices to the left (right) of the symbol (8> act on Dirac (taste) indices. 
We denote by 7 and t the matrices acting on these indices, respectively. The 
operators 

P^^^ = ^ [1 1 ± 7^.75 ® ht^.] . (71) 

are orthogonal projectors. The fermion fields are defined on blocks (see 
Appendix [O for details). The right and left derivatives V*^^) are given by 

V?) = ±1 (t(±) - 1) . (72) 
The factor 1/2 is due to the fact that the translate by one block 

rw) =5,,,,,±2A- (73) 

The model has a discrete chiral symmetry at m = 0: 

V'^-75®1V', -075(8)1. (74) 

To have an action bilinear in the fermion fields we introduce auxiliary scalar 
field ct(x), whose integration generates the four fermion coupling: 

^5' = E ' E '^(^) (rn + a + + ^ ^^(x) . (75) 

X y X 

The partition function now reads 

Z = j [(if7][##] exp [-S'] . (76) 

5.1 The effective boson action 

In the model with 4-fermion interaction with Kogut-Susskind fermions in 
the flavour basis [19] the matrix M = 0, while the matrix N is Hermitian 
and is given by 



Nia) = -2 < 



3 

(m + cr) 70 O 1 + E 7o7i (8 1 



p(-)vW+pWv(-)' 



(77) 
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We will restrict ourselves to flavour singlet composites, what means that 
^xK acts trivially on flavour indices, and thus, obviously, only flavour singlet 
composites can be written as linear combination of the $xK- From now on, 
we will ignore flavour indices. According to the results of the above Sections, 
the effective action at the saddle point is 



2^tr_ln y/l + H^ + H 
tr In 



t 

' 2 



VT+IP + H 



(78) 
(79) 



because the factor 1/2 coming from the projector Pq present in tr_, is 
compensated by the exponent 2 in the In. 

We now show that this is equal to the standard result. The latter is 
obtained by a direct integration over the fermion flelds in the partition 
function 



da] exp 



-Sp 



AN 



f 



2g' 



(80) 



where 



Sf 



-Tr ln(m + cr + Q) 
du; 



■tr 



tt/2 



2tt 



In 



+ -(1 
2 ^ 



cos 2uj) 



(81) 
(82) 



Here Tr (notice the capital case T) is the trace on all the entries of Q, 
namely time is included, and since Q acts on the whole fermion fleld there 
is no projector Pq \ The explicit integration over uj reproduces, apart from 
a constant term, the expression in ()79p . Therefore 



Seff(f,^» 



d_ 



S p 



Lr 



■tr 



a 



2 HVl + 



(83) 



5.2 Mass of the effective boson 

Let us consider the Gaussian fluctuations around the minimum of the action. 
To this end, let us write 



iK{x) 
a{x) 



Tk + ^k{x) , 
a + r/(x) . 



(84) 
(85) 
(86) 
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We can assume the real, because the are defined up to a global phase 
which can be eliminated by a redefinition of the ^xE"- 

Let us expand the eff'ective action Ses{i + ^* + ^-.'^ + fl) to second 
order in powers of (p* , and r]. The linear term vanishes due to the gap 
equation and the terms involving 93 and ip*'^* vanish on the saddle point. 
Ignoring the constant term S'efr(^ j^jO"), the quadratic part of the effective 
action reads 

S-eff = J]??(X)Z)(''^)(X,X')??(X') 

xK,x'K' 
x,x'K 

+ pUx)Df'^\x,x')r^{x') (87) 



where the D's are given by traces involving the matrices ^xXj ^xXj and N 
(see appendix [D] for details). Recall that the entries of these matrices are 
the fermion spatial coordinates and Dirac-taste indices, and the label K of 
the composite bosons includes Lorentz, taste and spatial (radial and angular 
momentum) quantum numbers. We expand the structure matrices <I?x_ft' in 
a basis of matrices H(x, n) which act only on the spatial coordinates of the 
fermions and which act only on the Dirac-taste indices 

^xA' = Y CxA:,xin7 S(xi, n) . (88) 

xin7 

We accordingly expand the boson fields 

iA:,xn7 V'ii'(i, xi) . (89) 

XI A' 

We consider two options for the matrices H(x, n) 

[^(x,n)]^^_^^ = -Jxi.x+in^xa.x-in' (90) 
[H(x,n)]^^_^^ = '5xi,x'5x2,n- (91) 

In the first case x has the physical meaning of the centre of mass coordinate 
of the composites and n of the relative distance of the fermions. Then the 
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matrix H must be separated into two contributions: one for x = 2rx,n = 
4r„, the other one for x = 2(rx + 1), n = 4r„ + 2, writing x = 2r, where the 
r are vectors with integer components. This accounts for the fact that the 
coordinate of the center of mass can faU in the fuh lattice, while the positions 
of the fermions and their relative distance are coordinates of blocks. In the 
second case x and x are the coordinates of the fermions. 
For the matrices T we require 



Tr {T\ 



5. 



77 



Note that the F^'s verify the identity 



p(~) p p(+) p 

-'^0 7 -'^0 — i 7 • 



(92) 



(93) 



If we are interested in the dynamics of the (/9-fields we must determine the 
structure matrices and therefore the coefficients c. But in this paper we 
have the limited goal of reproducing the spectrum of the bosonic composite 
r/. For this purpose we integrate out the (/9-fields and we do not need to 
determine the coefficients c. 

Integration of the (/j-fields generates an effective action for r/ of the form 



(x, x') rj{x') . 



(94) 

The explicit expression of the kernels appearing in this expression can be 
found in Appendix [Dl but we remark that they don't depend from the 
structure functions of the (/9-fields. Close to the continuum limit (ct — > 0), 
its behaviour is given by 



Ses[r]] 



^-^(-a;,-p)r?(u;,p) 



tt/2 J_7r/2 



(95) 



As o" — > 0, the integrals defining Zi^{a'^), Zp{a^), and W{a'^) as discussed in 
Appendix ini are dominated by the logarithmic infrared divergence, and we 
get 



Z^{a^) ~ Zp{a-') 



Ina^ + 0(1). 



(96) 



'~ 4 ~ (27r)2 

The logarithmic divergences are the same which appear in the exact solution 
of the model at large Nf, leading to logarithmic corrections to scaling and 
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thus to the triviahty of the continuum Hmit. In conclusion, the mass of the 
scalar is 



which is the exact relation with the chiral condensate obtained from the 
gap equation. Notice also that Euclidean invariance is recovered in the 
continuum limit: 



Summarizing, in this Section we have shown that the hypothesis of boson 
dominance and the associated formalism, applied to the 4-fermion interac- 
tion model reproduces the exact solution in the boson sector, namely, the 
correct gap equation, mass gap, and logarithmic corrections to scaling. 

We know that the spectrum of the model contains also one fermion of 
mass = a. This property is reproduced by an extension of the present 
formalism [16]. 

6 Summary and outlook 

In this paper we applied a general method of bosonization to relativistic field 
theories whose low energy excitations are dominated by bosonic modes. This 
is always the case in the presence of spontaneous breaking of continuous 
symmetries. One of our goals is indeed the study of low energy hadron 
dynamics in QCD. 

Under the condition of boson dominance the fermion partition function 
can be restricted to boson composites. To realize this restriction we use the 
formalism of the transfer matrix, which is close to the Hamiltonian formalism 
of nonrelativistic theories, and therefore very natural dealing with real or 
virtual bound states. The projection onto the subspace of composites in 
the partition function is realized introducing coherent states of composites. 
The projection operator is not exact, but its approximation is the more 
accurate the higher the index of nilpotency of the composites (the number 
of fermionic states which define their structure functions). Indeed it leads 
to exact results in the model with quartic interaction in the Nf oo limit, 
in which case the index of nilpotency also tends to infinity. Evaluation 
of traces in the fermion Fock space can be performed without any further 
approximation for theories which are quadratic in the fermion fields or can 
be put in such form, a condition satisfied by all renormalizable models in 
3+1 dimensions. It generates a functional form of the partition function 
with an effective action of the composites. In this derivation time and space 



m„ = 2a 



(97) 




1. 



(98) 



20 



are not treated in the same way, and Euclidean invariance must be checked a 
posteriori, but all other symmetries are respected including of course gauge 
invariance. 

We checked our theory on the 4-fermion interaction model in 3+1 di- 
mensions. Not only did we reproduce all the known results in the boson 
sector, but we have also been able to determine the structure functions of 
the condensed composite. It turns out that the polar factor of the struc- 
ture function is identical to that of the Cooper pairs of the BCS model of 
sup er conductivity. 

The study of the 4-fermion interaction model in our formalism is not 
complete, because we did not investigate if, in addition to the sigma, there 
are other bosonic modes. Indeed what we actually did was to prove that pro- 
jecting onto the composites subspace and then integrating over the fermionic 
fields gives the same results as a direct integration (without projection) of 
these fields as far as spontaneous breaking of chiral invariance and mass 
of the composite are concerned. We must also observe that while the gap 
equation does not depend on the number of space-time dimensions, our eval- 
uation of the mass of the sigma is restricted to 3+1 dimensions because we 
used completeness relations valid in this space. 

The effective action we derived can be used in the study of QCD by 
numerical simulations and, we hope, by a perturbative expansion along the 
lines of [11] and by keeping into account the subtleties related to the 1/-/V/ 
expansion [20]. In this case the effective bosons carry obviously the quantum 
numbers of the chiral mesons and chiral symmetry is broken by the conden- 
sation of the sigma-meson. The introduction since the beginning of the 
expectation value of the sigma-field as a variational parameter, should help 
in a numerical simulation, as advocated in Ref. [8]. But we must remember 
that the form factors of the composites depend in general on the fields of 
the elementary bosons coupled to the fermions, as we have explicitly seen in 
the model with 4-fermion interaction. For an actual numerical simulation in 
QCD therefore, we need a trial expression of the meson form factors, which 
should depend on the gauge fields, as also required by gauge-invariance of 
the effective action. Apart from the interest which it has per se, a pertur- 
bative approach might prove useful also to provide trial expressions of the 
form factors. 

We conclude by observing that the present formulation can be immedi- 
ately used also for studies of systems at finite temperature, while for finite 
fermion density we must first include states with nonvanishing fermion num- 
ber. 
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A Some notations and conventions 

We have used the following identities valid for Berezin integrals 




(101) 



(100) 



(99) 



If I a) is a fermionic coherent state 



a) = e 



|0) 



(102) 



then 



{a\a) = e"*"(0|0) = e' 



(103) 



and 




(104) 



Remark that 



u\a) = a \a) 



(105) 



which implies the relations 



(a/3|e-^«|75) 



(106) 
(107) 
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With the help of these formulae we will compute 

(a/3|e*^"e"^^^^^|0) (108) 



dYd-fd6*dS 



{aP\e^^^\jS){jS\e^^^^*^\0) (109) 



= J [dj*d-fd6*d6]e~^*'^-^'^+^^^+'''"^+^*^+'^''^^^' (110) 
= J [dS*d6] e-s^i^+^*^^)^+P*s-s*^*-* (111) 
= det+(l + 5*iV^)e"^*TT^^*"* (112) 

B The operator V 

The restriction of the partition function to the subspace of composite bosons 
can be written 

Zc = TT^l[VTt,t+i (113) 
t 

where V is the projection operator on this subspace. For the sake of sim- 
plicity we consider the case of a unique composite. In such a case 

7^ = ^— I (6t)"|0)(0|l>'*| (114) 

where 

i/„ = (0|$"|(4>t)"|0). (115) 

We must show that 

(0||.-7^(|.t)n|0) ~ (0|l>-|(|.t)"|0) = 6m,nl^m. (116) 

These equations are generated by the following ones 

imn ^ (e'lo (117) 

whose right and left hand sides are 

(^'i^ie") = /^exp(-£:(r,e,e'*,e")) 

i^'ia = exptr+ln(l + e'*?"$$t) (118) 
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where 



C C") = tr+ ln(l + CC^^^) - ln(l + C'*?^^^) - ln(l + C^"^'^^) 

(119 

We evaluate the integral by the saddle point method. The saddle point 
equations are 



(r-e'*)tr- 



with solutions 

At the saddle point 



(1 +^*^$$t)(l +^*^"$$t) 
' (1 +^*^$$t)(l + ^'*^$$t) 
^ = ^" ^* = 



^^r , e, r , n = -tr+ in(i + c'*e"<f $ 



Moreover 



In conclusion 



tr. 



;i +^'*C"$$t)2 



tr_ 



(120) 
(121) 
(122) 

(123) 
(124) 

(125) 



In presence of a number of composites n << the desired result, together 
with the idempotency property of projector 



follow if we assume 



-p ~ "P 



tr l^^'^'Y ~ 



-n+l 



(126) 



(127) 



But in states with n ~ fi, it is not possible to satisfy (jll7p even with 
an absolute freedom about the form of the structure functions (which are 
instead determined by the dynamics). The best we can do is to satisfy (jllTp 
(apart from an irrelevant constant factor) for states with n + k composites, 
for fixed n ~ and |A;| << il. 
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To clarify this point let us go back to the case of many composites and 
consider their commutation relations 



V. 



(128) 



In states with a number of composites n « fl, these relations are approx- 
imately canonical provided the structure functions are sufficiently smooth. 
Indeed in such a case the last 2 terms are of order n/Q,. 



C The transfer matrix for Kogut-Susskind fermions 

Kogut-Susskind fermions in the flavour basis are defined on hypercubes of 
twice the ordinary lattice spacing. We shall have a fermionic field "(t,x) 
(and of course the corresponding ip), where i = {1, . . . ,Nf} is the fiavour 
index, a = {1,...,4} is the spinorial index, a = {1, ...,4} is the taste 
index, while (t, x) is a 4-vector of even integer components ranging in the 
intervals [0, Lt — 1] for the time component while [0, Lg — 1] for each spatial 
components. 

In order to simplify the notations let restrict here to periodic functions 
defined on one dimensional interval x G [0, L — 1] , the extension to many 

dimensions and anti-periodic functions being obvious. 

The sum over even sites x includes for convenience a factor 2 

E'-2^- (129) 

X X 

Momenta are quantized according to 

2vr L , , 

p=—n n = 0,l,--l (130) 

The Fourier series and its inverse for the function f{x) defined on even sites 
are 

L-2 

hp) ■■= 5^7We^"^ (131) 

x=0 

L/2-1 , , 
n=0 ^ ^ 

Notice that we use a double tilde to distinguish this transform from the 
standard one for a function defined on all the sites which will be denoted by 
a single tilde. As / is a function of an even argument we could set 



g{s) := /(2s) (133) 
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where s is an integer in the interval [0, L/2 — 1], then their Fourier transforms 
are simply related as 

L/2-1 

hp) = 2 ^ /(2.)e^2.p 

s=0 
L/2-1 

= 2 9{s)e'^'P (135) 



2g{2p) (136) 



In the infinite volume limit we get the Fourier transform 

oo 

hp) ■■= (137) 

fix) := / //(P)e-'^^ (138) 

As a consequence, for instance, the 5-function must be written 

6{x -y) = 2 r''^ ^ e*P(^-^) (139) 
J-n/2 2vr 

and the Fourier transform of a matrix of the form 

A = F(-A) (140) 

for an arbitrary function F is 

A(p) = 2F Q(l - cos 2p)^ (141) 

and its trace 

(142) 

We report here the expression of the transfer matrix in the flavour ba- 
sis |19] . because the expression in the spin-diagonal basis is not known to 
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us in a convenient form. The matrix N, which is Hermitean, in presence of 
a gauge field denoted by C/, is given by 

N{a, U) = (143) 



-2 ^ (m + cj)(7o 1)+ Y^ilolj ® a) 

i=i 



p(-)vW + pWv(r) 



where 



Vf = 4(f/f^)(t)r,W-l) (144) 



are right-left covariant derivatives. The [/(t)'s are operators whose matrix 
elements are the Wilson link variables ?7^(t,xi) 

^ Xi,X2 

The matrix M depends on the link Uq variable. In the gauge Uq ~ 1, in 
which Uq = 1 with the exception of one time slice, M = 0. If we adopt this 
gauge, we must impose the Gauss constraint on the states. But it is very 
simple to write the expression of M without gauge fixing following Ref . ^21j . 
This is actually necessary in a numerical simulation. 

In the absence of gauge fields and with constant a = a, according to 

m 

= {m + af - A. (146) 

with 

^ = I E + - 2) (147) 

1=1,3 

The eigenvalues of are therefore the fermion energies 



where the momentum q is 



E^ = {m + of + , (148) 



3 



t = Y.^^ (149) 



i=l 



with 



g2 = i(l_cos2(7i). (150) 
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D Quadratic fluctuations 

Let us give some details about the expansion of the effective action around 
the constant fields {B,B\a) which minimize it. We shall set 

Bt = B + 5Bt (151) 
C7(t,x) = a + r/(t,x) (152) 

where 



B = ^* • $ (153) 
6Bt = i^'' ■^)t = Y,VK{t,^T^^K (154) 



and shall concentrate on quadratic fluctuations. 

In the parametrization of the form-factors <I>xX we choose the form 
in (j9ip for the matrix H(x,x'). 

It is also convenient to introduce the matrix fjt by 

r?t = 2 ^ r/(t, x) H(x, x) (70 ® 1) . (155) 



D.l Expansion of the action 

According to our definitions the expansion for the effective action reads: 

S,s{v\^, r,) = 5eff + tr+ { R (B6bI + 6BtB^) 

t 

- Rn{Bn5B\^^ + 5Bt B^N + BnvI+2 + VtB^ 

+ R 6Bt 6bI - Rn {SBt 8B\^^ + fitfjl^^ + rjt SB^^ + 6Bt 4^^) 

-^[R(B6B^ + 6BtB^)f 

+ ^[Rn (Bn 6b1^2 + % + BNfil+2 + m%) f } , (156) 

Bn = N + B (157) 
R = (^1 + SS^) ^ (158) 

Rn = (i + BnB^nY^ ■ (159) 



where 
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Ignoring the constant term, the expansion can be cast in the form 

+ ^ <^|, {t, x) ipK' {f, x') Dfj,^^ {t, x; i', x') 
txK,t'x'K' 

+ Yl »?(t,x)¥.i,(t',xOZ)i^'^)(t,x;t',x') 
tx,t'x'K 

+ Yl v{t,^)<PKit'y)D^K''''\w'y), (160) 



tx,t'x'K 



since the linear terms vanish due to all stationarity equations — we are ex- 
panding around a minimum — and the terms involving '^k'-PK' and 
vanish due to the stationarity equations for the fields B and . 
The D's entering the above equations are: 

L>(^'')(t,x;i',x') = 4tr+i?jv {-2(5t.,t+2H(x,x)(7o0l)H(x',x')(7o® 1) 

+ 5t',t+2 S(x, x)(7o (8) 1) B^NRNBN^i^', xO(7o «> 1) 

+ St^t'+2 BjvH(x,x)(7o 1) i?ivH(x',x')(7o 1)-Bjv 
+ (5tt/ [BnE{x, x)(7o «) 1) RnBnE{:>c', x')(7o «> 1) 

+ S(x,x)(7o®l)Bjvi?ivS(x',x')(7o®l)sjv]} (161) 

L>S^;.^^ (t, x; i', x') = tr+ [Sw R^^k [l - B" Rb) 

-5t>,t+2 Rn^xK (l - B^RnB) (162) 



D('''^)(t,x;t',i^) = 2tr+i?^ {-,5t.,t+2S(x,x)(7o®l)4^ 
+ 5tt' S7v$ix-RjvS;vS(x, x) (70 ® 1) 
+ <5t,,t+2S;v$iKS(x,x)(7o®l)5jv} (163) 



D^^^*\t,x;t',K) = 2tr+i2^ {-5t,t,+2^xxS(x,x)(7o®l) 
+ (5tt/ $xifS!v^JvS(x, x) (70 (8) 1) 

$xi^5]vi?ivSivS(x,x)(7o 1)}. (164) 



+ 
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D.2 Evaluation of the traces 

We need to make explicit the meaning of the indices K^s : K = (11,7), so 
that 

<^K > ¥'n7- (165) 

and therefore 

^ (166) 

Now we can readily compute the traces entering equations (|16ip - (|164p . 
In order to present the results, it is convenient to introduce additional defi- 
nitions. 

On the solution of the gap equation, the matrices R and Rm can be 
written as 

R = 4^ (167) 
Rn = (168) 

where 

A = . (169) 

Obviously, R + Rjq = I. Notice that the following relations are verified: 

= (4i?)-2 _ A, (170) 
{l + Af = {4H)-^ + 1 + A. (171) 

We also introduce the following notation: 

e = RNil + A), (172) 
Sf) = G(i;(')-1), s = +,-, (173) 



and 



r = tr'_[(7o»l)r^ ], (174) 
tl = tr'_[(7,7o0ll)P/-^)r^], (175) 
= tr'_[(7,®Il)pj~^')(7o7i®l)J^^"'^r^ ] (176) 
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where tr'_ is the trace on the Dirac-taste space. Introducing the Fourier 
transforms 

^n,«.x) = /^/^/|^e'"e*-e'<.-'^„Kp) (178) 
we get the following effective action: 

^"=1"^! i^5^<'"(-,p)^(--,-p)^(-,p) 
^IyJw?! ||^r'(-,P)^i,(-,P),(-.P), (ITS) 

where 

i^('"')(a;,p) = 1 + 64y"^ |-(2-cos2u;)i?Ar(q)i?iv(p-q) 
- ^ E[^i^^(q)^i"^(P - + s\-\d)s\'-\v - q)] 
+4a'e(q)e(p + q)| , (180) 
^j;^ (^,P) = '^77' [^(P) ^(q) - e*'"^iv(p) ^iv(q)] , (181) 

L>fe^^(c^,p) = 2 [4a2e(p) e(q) - e^2-it;,v(p) i?jv(q)] 

- 2a ^ .ij, [G(q) 5f)(p) + e(p) 5l-^)(q)^ 
j,s=±l 

+ E E -'C.v5^^q)^r(p)|- (182) 

i,i s,s'=±l I 
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D.3 Effective action for the //-field 

The integrations over ip^^ and </?*^ give the following effective action for 77: 

where 

A('"')(a;,p) = 



(183) 



q 



D. 



{■<rf) ' 



(--^S P - q) -Dqf ■ (^', p - q) 



(27r)3 i? (q) i? (p - q) - e^Su; (q) (p - q) ' 

(184) 

The sum over 7 appearing in the above equations can be performed using 
the following identities 



7 

.7*,7 

7 

E/7*+7 
''is ''js' 

7 

E,7* ,7 



tr'_ [loll 



tr'_ 



0, 



= tr'_ 



tr'_ 



(7, ®l)Pf ^(7,®l)^•^■ 



4 5ij Sss' 



tr'_ 
0, 



(7,- C5 1) ^ (7fe C5 1) ^ (7i ® 1) p. 



(185) 
(186) 

(187) 



(188) 

i-s) ' 

(189) 



A, 



.7 



= tr'_ 



(7., a) p/f ^ (7,2 ® a) 4''^ ^1 ''^ (^^1 ^ 1) ® 1) 1 (190) 



which follow from the completeness relations satisfied by the matrices F-y's 
which form a basis in the relevant subspace of matrices acting on Dirac and 
taste indices. 
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With the above equahties it is easy to get the exphcit form of A^^''^ ('^, p): 
A('"')(a;,p) = (191) 
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f rf^q 

J (27r)3 i?(q) 



) (p - q) - Rn (q) Rn (p - q) 

X I [4^2 e (q) G (p - q) - e'^'^RN (q) i?iv (p - q) ]' 

- [4a^e (q) e (p - q) - e^^'- Rn (q) i?7V (p - q) ] 
X ^ [ (q) (p - q) + St^ (q) sj") (p - q) ^ 

i 

+ ^ (q) (q) ^ (p - q) S^^ (p - q) 

j i 

+ 4 ^ I e (p - q) 5f ) (q) + G (q) (p - q) | ' | 

Now we know D^^^\u),p) and A^^^i\ui,p) in terms of integrals of known 
functions of momenta and a. 

D.4 Low energy expansion 

Let us consider the expansion of D^*"?) (a;, p) and A(w) (o;, p) in powers of lo 
and p up to second order: 

L>(''^)(u;,p) + A('"')(a;,p) ~ Z^{w'')iu^ + Zp(w'')p'' + ^(a^)^^ (192) 

The gap equation impUes that the above expression at = and p = 
vanishes close to the critical point, that is as a — > 0. It is straightforward 
to get the term at a; = p = 

W(a') = 16 / ' ^ a' (193) 

that is 

At a; = 0, by keeping the most singular contribution when ct — ^ 0, we get 
the contribution proportional to p^ 



2 da 1 
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Similarly, expanding in u; ,at p = 0, since the linear term is not diverging, 
the dominant contribution is proportional to iJ^ and 
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